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The notions of pure states and inherent structures, i.e. stable configurations against f-spin flip 
are discussed. We explain why these different concepts accidentally coincide in mean-field models 
with infinite connectivity and present an exactly solvable one dimensional model where they do 
■ not. At zero temperature pure states are to some extent related to fc-spin flip stable configurations 

with k — > 00 after the thermodynamical limit has been taken. This relationship is supported by an 
explicit analysis of the TAP equations and calculation of the number of pure states and fc-spin flips 
stable configurations in a mean-field model with finite couplings. Finally we discuss the relevance 
of the concepts of pure states and inherent structures in finite dimensional glassy systems. 
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I. INTRODUCTION 



S3 . 

In a series of seminal papers, Kirkpatrick, Thirumalai and Wolynes (KTW) suggested ten years ago that (discon- 
tinuous) mean-field spin glasses could serve as a paradigm for vitreous systems The major success of KTW theory 
lies in its ability to predict the existence of both low (Kauzmann) and high (dynamical freezing) temperature critical 
points and thus to describe in a unified framework the thermodynamical and dynamical signatures of the glassy 
transition ||. This achievement relies on a simple but rich picture of glassy pure states (PS), that is local minima in 
the free-energy landscape, and of the dynamical evolution taking place in the latter. To what extent KTW mean-field 
1 . picture applies to real systems is a crucial issue . 

An apparently related approach based on the investigation of the potential energy landscape and of its local minima 
called inherent structures (IS) had been already proposed at the beginning of the eighties by Stillinger and Weber 
in the context of liquid theory ||. IS present a considerable advantage with respect to PS: they are well defined 
and free of any mean-field hallmark. Recently, numerical studies proposing that some understanding of the glass 
\ transition could be gained from the analysis of potential energy landscapes Q strengthened the feeling that Stillinger 
^ . and Weber's analysis and KTW theory were basically the same description of glassy systems. This point of view was 
later supported by the equivalence of IS and PS in some mean-field spin-glasses, e.g. Sherrington-Kirkpatrick (SK) 
or p-spins models and by analytical works |||| . 

In this letter, we present some arguments to clarify the relationship between IS and PS. We present some one 
I . imensional models for which IS can be exactly computed and do not coincide with PS. We show in addition that 
Q\ ' the degeneracy between IS and PS in infinitely-connected models is due to the vanishing of the interactions in the 
thermodynamical limit and is lifted by the introduction of finite couplings. Focusing on the TAP equations "defining" 
PS and the number of the latters, we show that zero temperature PS are related to extended IS, that is configurations 
stable with respect to an arbitrary large number k of spin flips (and not only to a single one as usual IS) §. 



II. DEFINITIONS AND NOTATIONS 



a 

o 

We consider a model including N Ising spins Si . A configuration of spins will be said fc-stable if its energy cannot 
be decreased by flipping any subset of k (or less than k) spins. Let us call s&(e) the logarithm (divided by N) of the 
number of fc-stable configurations with energy density excess e — ecs with respect to the ground state energy density 
eos when N — > 00. IS correspond to 1-stable configurations and si(e) is usually called configurational entropy. We 
' have of course Sk+i(e) < Sfe(e) and Sfc(e) = —00 for e < eos- 

The notion of PS is much trickier. In mean-field systems PS are usually defined through TAP equations [[l0| for 
the local magnetisations m^. TAP solutions are exponentially numerous in the volume N and can be accounted for 
through the so-called complexity er(/), that is the normalised logarithm of PS having free-energy densities /. We 
shall focus in the following on the zero temperature limit of the complexity <r(e) and compare it to Sfc(e). 
We shall denote hereafter the maximum over e of s/c(e) (respectively a(e)) by Sk (resp. a). 
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III. ONE DIMENSIONAL DISORDERED CHAIN 



We first consider a Ising chain with disordered nearest neighbour couplings Li and no external field The LiS 
are independent random variables drawn from an even distribution P(L). No phase transition can take place at finite 
temperature: there is only one paramagnetic PS 111]] . Therefore, the zero temperature limit of the complexity simply 
reads a(e) = if e = eas, = — oo if e 7^ e G s where the ground state energy density equals ecs — — j dLP(L)\L\. 

Ettelaie and Moore have calculated si(e) [Q. We extend their approach to obtain s&(e) using the following 
observation by Li: a necessary and sufficient condition for a configuration to be fc-stable is that each frustrated bond 
be fc-weak, i.e. smaller in magnitude than the 2k nearest bonds (k on the left, k on the right) on the chain The 
typical and self- averaging Sfe(e) thus reads 



Sfe(e) 



lim — In 




e-as 




(1) 



where Nf. is the number of fc-weak bonds on the chain and Tg equals 1 if the I th fc-weak bond is frustrated, otherwise. 
Using an integral representation of the delta function in (Q), the r/s may be traced over and one finds 



^ r p+ioo 
Sk (e) = lim — In < / 
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2nNi 
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where Pk{L) = P(L)(2 J^dL' P(L') ) 2k denotes the probability that a bond is fc-weak and has value L. The 
integral over u can be evaluated by the method of steepest descent. In fig.l, s/-(e) is plotted for a Gaussian P(L) 
(with variance unity) and for different values of k. At fixed k, Sfe(e) is symmetric around its maximum located in 
efe = ggs + / dLPk(L)\L\, Sk = ln2/(2fc + 1). Note that the height of the maximum Sk, related to the total number 
of fc-stable configurations is independent of P as found by Li |Is)| . On the left boundary, all curves exhibit an infinite 
slope when e — > ecs- More precisely, Sk{&GS + e) oc — elne if the support of P(L) does not contain L = and 
s k {e G s + e) oc e («+ 1 )/( Q + 2 ) if p(L) oc L a when L -> (with a > -2 to make e GS finite). 



s k (e) 




FIG. 1. Logarithm of the number of fc-stable configurations Sfc(e) as a function of the excess of energy density e - 
Gaussian coupling distribution P(L) (with variance unity and zero mean). From top to bottom: fc = 1, 2, 3, 4, 5. 



bgs for a 



The above model is a clear example of a system with many IS and at the same time a single PS. More precisely, 
the zero temperature limit of the complexity er(e) (related to PS) is not equal to the configurational entropy s\(e) 
counting IS but to Soo(e). Before exposing why s\ and a accidentally coincide in infinitely connected spin-glasses, let 
us briefly see to what extent the above results apply to generic finite dimensional systems. 



IV. GENERAL REMARKS FOR FINITE DIMENSIONAL SYSTEMS 



The main features of Sfe(e) can be found by some general arguments which apply to any finite dimensional system. 
The behaviour of si(ecs + e) at small e has been predicted in JT^JT^] and depends only on the lowest-lying excitations 
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of the system. In addition, the scaling of the right-edge energy 2ek — ecs at which Sfe(e) vanishes can be understood 
by a simple energy balance argument. Flipping a cluster of k spins among a configuration with energy density e, 
one may gain a bulk energy Eb(k) ~ — k(e — eos) and loose at least some surface energy E s (k). As a consequence, 
fc-stable configurations with energy densities higher than ecs + E s (k)/k should not exist. For 1-D Ising disordered 
chain, E s (k) is equal to the average magnitude of the smallest bond among k ones and can be easily computed. As 
expected E s {k)/k ~ — eos for large k. Let us emphasise that the previous argument implies that there cannot 
be oo-stable configurations (after the thermodynamical limit has been taken) with energy densities differing from the 
ground state one: Soo(e ^ ecs) = — °°- 



V. MEAN-FIELD SPIN-GLASSES WITH INFINITE CONNECTIVITY 

It is well known that the zero temperature limit of the complexity equals the configurational entropy in infinite 
connectivity spin-glasses |l|] . For instance in the case of the SK model, 

a SK = sf K = max j-y + In 2 J Dy\ ~ 0.1992 (3) 

with Dy = e~ y2 1 2 / \f2n dy |l^Jl^]. The equivalence between TAP solutions at zero temperature and 1-stable configu- 
rations (IS) is a straightforward consequence of the vanishing of couplings in the thermodynamical limit. The scaling 
of the couplings Jy with size N ensures that the excitation energy AE iy corresponding to flipping spin i remains 
finite in the thermodynamical limit, e.g. Jij = 0(l/yN) for the SK model. As a consequence the excitation energy 
AEij = AEi + AEj — 2JijSiSj corresponding to flipping both spins i and j reduces to AEi + AEj in the large N 
limit. For infinite connectivity models IS are also fc-stable configurations for every finite k. Therefore, Sfc(e) does not 
depend on k and s± = . . . = s^, = a. 



VI. TAP SOLUTIONS AT ZERO TEMPERATURE ARE oo-STABLE CONFIGURATIONS 

We now show that in presence of finite couplings, the degeneracy between all s^s is lifted since A;-stable configurations 
are not necessarily k + 1-stable. In addition, we show that zero temperature TAP solution are oo-stable and not only 
1-stable. To do so, let us first recall how TAP equations can be derived for a one dimensional spin glass [^9|. We 
call U and respectively the effective magnetic fields due to all the spins on the left and on the right of i and hi the 
external magnetic field acting on Si- It is easy to write the equation verified by rrii in terms of U-\ and r^—i, 

rrii = tanh [(3h n + tanh _1 (tanh(/3Zi_i) tanh(/3Jj_i)) + tanh _1 (tanh(/3r l+ i) tanh(/3Jj))] . (4) 

This is the equation verified by the local magnetisation of a spin Si interacting with two spins <Si_i, Si+i on which 
the magnetic fields and r,_(-i act. Focusing on two neighbouring spins the equations verified by the two local 
magnetisations are: 

tanh(/%) + tanh(/3Jj) tanh(/3r,-+i) . 

(5) 



1 + tanh(/%) tanh(/3Ji) tanh(/3r j+ i) 
tanh(/3r i+1 ) + tanh(/3J i ) tanh(/3ij) 
mi+1 ~ 1 + tanh(/3^) tanh(/3 J t ) tanh(/3r l+ i) ^ ' 

These equations give m, and mj+i as a function of U and r^ + i. Inverting (^), (^|) one obtains U and r^ + i as functions 
of rrii and m^i. Plugging = ij_i(mi_i, mi) and r^i = r,+i (ra;, TOj+i) into (^) the TAP equations for the local 
magnetisations are established [^9) . 

To prove that TAP solutions are mapped into ground states in the zero temperature limit let us focus on k 
contiguous spins on the chain. The TAP equations on these k spins are by construction the equations verified by the 
local magnetisations of the system of k spins, as if it were isolated from the rest of the chain and there were magnetic 
fields li_ , ri + acting on the leftmost and rightmost spins respectively. As a consequence in the zero temperature limit, 
the k local magnetisations tend towards the configuration of k spins which realizes the global minimum of the isolated 
system of k spins for any value of U_ and n + . As k can be made arbitrarily large once N has been sent to infinity, 
the configurations corresponding to the zero temperature limit of TAP solutions are oo-stable. The argument can be 
straightforwardly extended to Ising spin glasses on Cayley trees or on random graphs. 
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VII. CALCULATION FOR A MEAN-FIELD MODEL WITH FINITE COUPLINGS 



The results of the previous sections suggest that in finite connectivity models Sk is a non trivial function of k which 
tends towards the zero temperature limit of the complexity a for large values of k. To corroborate this point, it 
is interesting to study a model whose complexity a is non zero. Consider the Hamiltonian H = — ^2 i< j JijSiSj — 

12iLi LiSiSi+i . The first term is the usual SK Hamiltonian: couplings are independent Gaussian random variables 
of zero mean and variance J 2 /N. The one dimensional Hamiltonian involves disordered interactions L$ which are 
independent random variables vanishing with probability 1 — a and distributed according to an even distribution law 
P{L) with probability a. This system smoothly interpolates between the SK model (a = 0) and the one dimensional 
spin glass (J = 0). We calculate below the departures of a and Sk from their common value in the SK model to the 
first order in a only to avoid useless tedious calculations. 

To compute the complexity two approaches can be used. The first one consists in making the (annealed) sum over 
all locally stable TAP solutions ]l7| ]. The other method relies on the computation of the logarithm of the (average) 
n th power of the partition function using Bray and Moore replica symmetry breaking scheme pp| ]. As n goes to zero, 
the annealed complexity is recovered pl| . Calculations are sketched in Appendix A and B and provide the same 
average complexity a(a) = a SK + a 0(a 2 ) where 

f+00 />+oo 

„(1) ' ' 




7 W = -1+ / dLP(L) Dy 1 Dy 2 9(y 1 +y 2 + 2z*)} / { I Dy ) . (7) 

\J — 00 J — z*—L/J 

and z* ~ 0.506 is the maximum of (||). 6 denotes the Heaviside function. 

We now turn to the computation of Sk- The excitation energy corresponding to k spin flips depends on the SK 
couplings and the initial configuration through the variables Xi = JijSiSj only. For large N, the XiS become 

independent of the spin configuration {Si} and can be written as x, = yi + z, where j/j and z are independent Gaussian 
random variables of zero means and variances J 2 and J 2 /N respectively [£2| . The annealed configurational entropy 
is defined through 



e 



N 

ATsi _ 



= (J2I[H^iS^iS l +L l S t S t+1 +y l + z)) . (8) 



Once the average (•) over disorder has been carried out, spins may be traced over and (||) may be evaluated for large 
iV by the method of steepest descent. We find s\(a) = a SK + a + 0(a 2 ) with 




+ 00 r+00 \ / / P+OC 

dLP(L) / D Vl Dy 2 / Dy ) , (9) 

-00 J -z*-L/J J I \J -z* 



giving si > a as expected, compare (Q) and @ p3j . Repeating the same calculation for k(> 2)-stable configurations, 
we have found that they are equally numerous to the first order in a: Sfe(a) = a SK + a Cjp + 0(a 2 ) with c£ = 

in) 

for any k > 2. This result extends to higher orders in a as follows. Let us call c k the coefficient of a n in s/.. 

(n) 

Contributions to c k come from disordered lattices of L bonds involving clusters of at most n + 1 contiguous spins 
along the chain. Clearly, spin configurations on such lattices that are n + 1-stable are for ever stable. Therefore, 
> c 2™^ > • • ■ > c%+i = Cfc^ f° r an y k > n + l. We have in particular explicitely checked that > 
As a conclusion, we have found a mean-field model with non zero complexity a and si > a. We have checked that 
a coincides with Soo to the first order in the density of O(l) couplings. 



VIII. CONCLUSION 



The above examples have shown the differences between IS and PS even in mean-field models provided that finite 
couplings are present. Our results could also be extended to models with continuous degrees of freedom, e.g. systems 
of N interacting particles mJ. In this case IS are configurations stable against infinitesimal moves of any subset 
of the N particles, that is local minima of the potential energy, fc-stable local minima can be defined as IS stable 
against finite moves of any subset of k (or less than k) particles, oo-stable minima gather all configurations that can 
be reached from each other by crossing finite barriers and are related to PS. 
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From the equilibrium point of view, PS and thus complexity are the only relevant concepts. The decomposition 
of the partition function put forward by Stillinger based on IS though mathematically exact does not seem to be 
thermodynamically founded j^J. However, usual nucleation arguments imply that no PS with free-energy / ^ f eq 
can live forever, that is a(f ^ f eq ) = — oo. We thus face the following alternative. Either nucleation arguments break 
down for some reason and finite-dimensional glassy systems with full curves for a(f) may exist Or they hold 

and the KTW entropy crisis scenario cannot be extended to realistic systems without taking into account the notion 
of life time for PS (as was already emphasized in pi). IS, besides their phenomenological interest can be seen as a 
very valuable step in that direction. 

Acknowledgements: We are very grateful to Marc Mezard for very interesting and motivating discussions. 



IX. APPENDIX A 



To first order in a spin i has at most one neighbour v = i ± 1 along the chain of L bonds. As a consequence the 
TAP equations can be written as 

m < = l + uXnZ L L) ' ^^E^i-W-^ , (10) 

if? 4 *) 

where q = 1 / 'N X^iLi m i- To compute the annealed complexity we sum over all solutions of equations jTc|) following 
fl7f . Introducing auxiliary variables one can perform the average over and using the same notation as |l7| , we find 

a = extr 9)A ,A {~M - A(l - q) - A 2 /2(3 2 J 2 + (11) 
1, / f+'n dtj 1 / (tanh -1 U - Am,) 2 , 2 \ \ \ 



I {L}. 

Performing the average on Li at the first order in a and taking the zero temperature limit we obtain the expression 
of 7. 



X. APPENDIX B 



To compute the n th moment of the partition function, we use the replica trick and obtain n coupled chains. Defining 
the overlap q ab between replicas a and b and the 2" x 2" transfer matrix, 



M(S,T) = {l-a + a / dLP(L)exp 



f3L S a T a 



exp 



B 2 J 2 . 

\ „a& / nfl nb i rparnb\ 



a<b 



q ab (S a S b + T a T b ) 



(13) 



the complexity may be computed within Bray and Moore two-group Ansatz [j20|: a = extr{ g }[— (/3J) J2a<b(l a ) /2 + 
In A({g})] where A({q}) is the largest eigenvalue of M. The saddle-point equation for {q} read q ab = J2s S a S b [^(S)] 2 . 
The normalised ground state wave function ^ obeys the two group symmetry with breakpoint m l20j] . $ depends 



on 



S = (S\...,S n ) through S_ = T, a < m S a and S+ = E a >m 5,0 onl Y and can be written as 



c. J dhadhiplho, hi) exp(f3ho(S+ + S-) + hi(S+ — S-)) where c is a normalisation factor. In the zero temperature limit, 
the effective fields ho, hi are both of order one and the overlaps q ab read q± = A — z/((3mJ), q2 = A — w / ([3m J) 2 , q% = 
A + z/(/3mJ),A — ► 1 with the notations of 
effective field distribution 



The eigenvalue equation for A, provides a linear equation for the 



Xp(h , hi) 



dLP{L){ (1 - a)S(L) + aP{L) 



DuDv J dhgdh'i exp yh!iWL(h' Q 



>pyi — Ju, h'i ~ uz — v\J 2w — z 2 ^j S yho — </?l(/iq)^ 6 yhi — ^i^C^d! 



(14) 



where A = Aexp(/3Jmz). ul(x) (respectively lpl(x), tjl(x)) equals —1 (rcsp. — L, 0) if x < —\L\; (resp. x sign(L), 
sign(L)) if — \L\ < x < \L\ and 1 (resp. L, 0) if x > \L\. The complexity a is found through maximising — z 2 /2 — w + 



5 



In X(z, w) over z, w and can in principle be computed for any a. The small a expansion of a is computed from the 
expansion of p around the SK distribution p SK (ho, hi) — 8(ha)5(hi) using (|l4|). 
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